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1. Introduction 


Low-frequency airborne radar systems are a promising technology for the 
surveillance of foliage-obscured stationary and moving ground targets such as 
humans and vehicles. 1 ^ Although much effort has been put forth in this area over 
the years, many unsolved problems remain as the foliage ecosystem itself—with 
scatterers of various shapes and sizes in multifarious arrangements—is inherently 
highly cluttered, and therefore, poses an intricate challenge for target detection and 
recognition, especially for a sensing platform in motion. Determination of target 
detectability ultimately necessitates a thorough fundamental understanding of the 
target and foliage scattering phenomenology (as well as the tree clutter statistics 
and attenuation) as a function of frequency, polarization, signal incidence angle, 
and foliage properties. As such, a full-wave electromagnetic solver is proposed in 
this study for scattering and propagation characterization in the foliage 
environment. Specifically, the analysis framework employs a parallelized finite- 
difference time-domain (FDTD) algorithm for deriving the far-field responses of 
complex realistic tree structures at P-band (200-500 MHz): the backscattering 
behavior of a single tree is first considered, and then large-scale simulations are 
performed to examine the responses of a forest stand in both the frequency and 
imaging domains. The solver is intended as a computational tool for evaluating the 
effectiveness of detection modes based on methods such as airborne synthetic 
aperture radar and ground moving target indication; nevertheless, the results and 
observations featured are also expected to be useful for remote sensing studies in 
which the primary interest lies in retrieving forest parameters from the scattering 
return. 

2. Simulation Framework 


As the first step, 3-dimensional representations of realistic tree structures need to 
be systematically constructed. The geometrical models in this work are created 
using the open-source random tree generator Arbaro, 5 which implements the rule- 
based growth algorithm developed by Weber and Penn. 6 Essentially, each tree 
structure is generated by following a series of recursive rules derived from 
geometrical observables. The tree structure in its entirety (trunk, multi-level branch 
complex, and leaves) can be exported from Arbaro in mesh form. Given the long 
wavelength at P-band, the effects of the leaves are often assumed to be negligible; 
therefore, only the trunk and branches are retained here for the electromagnetic 
simulations. Shown in Figs. 1-3 are some tree model examples. 
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Fig. 1 Quaking aspen {Populus tremuloides) tree model with increasing structural fidelity 
from left to right 



Fig. 3 Eastern cottonwood (Populus deltoides ) tree model 

After the mesh has been properly processed (through operations such as closing of 
holes, translation and scaling, etc.), voxelization of the structure is performed with 
the parity count and ray stabbing methods. 7,8 The scattering responses of the 
resulting voxel grid are then calculated with a customized far-field FDTD 
algorithm. 9,10 The implemented code is fully parallelized and runs on large 
distributed computer systems using the message-passing interface (MPI) 
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framework. In effect, the simulation domain is decomposed into rectangular sub- 
domains and the FDTD equations are solved separately for each sub-domain within 
one MPI process. The computational approach is highly scalable, even when the 
number of MPI processes is in the hundreds. More details on the electromagnetic 
solver can be found in our papers. 9,10 

3. Numerical Experiments 

The backscattering return of a standalone tree above a half space is first analyzed. 
The tree species considered is the quaking aspen (.Populus tremuloides), as shown 
in Fig. 1. Both the tree and the ground are assumed to be electrically homogeneous, 
with relative dielectric constant and conductivity ( s r , c d ) of (13.9, 39 mS/m) and 

(5.45, 20 mS/m), respectively. The tree has a height of 7.4 m, with a base diameter 
of 34 cm and a canopy diameter of 3.4 m. The co-polarized responses (radar cross 
section and phase difference) derived from simulating the structure with different 
levels of geometrical detail are displayed in Figs. 4-6, as a function of frequency 
and elevation incidence angle Q t . The results are obtained by averaging the 

responses over 36 azimuth incidence angles <f ) j . For the set of parameters employed 

here, it is seen that the backscattering mechanism of the trunk-ground dihedral 
provides the dominant signal response. The branch canopy has a tendency to 
enhance the radar cross section; this is more evident and consistent at horizontal- 
horizontal (hh) than at vertical-vertical (vv); the inclusion of the secondary branches 
in this case, however, does not significantly affect the overall return. Also note that 
the impact of the branch canopy seems to decrease as the incidence angle 
approaches grazing. In general, the vv responses on average are lower—and more 
sensitive to the incidence angle—than the hh responses; both of these observations 
can be attributed to the Brewster angle effect. 
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(a) 



(b) 


Fig. 4 Average response (co-polarized radar cross section and phase difference) of a single 
tree as a function of frequency and structural fidelity at 0 t = 45° . Dotted line = trunk; dashed 
line = trunk + primary branches; solid line = trunk + primary and secondary branches. 
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(b) 

Fig. 5 Average response (co-polarized radar cross section and phase difference) of a single 
tree as a function of frequency and structural fidelity at 6 t = 60° . Dotted line = trunk; dashed 
line = trunk + primary branches; solid line = trunk + primary and secondary branches. 
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Fig. 6 Average response (co-polarized radar cross section and phase difference) of a single 
tree as a function of frequency and structural fidelity at 0 i . = 75°. Dotted line = trunk; dashed 
line = trunk + primary branches; solid line = trunk + primary and secondary branches. 

Next, the simulation of a large forest scene consisting of 36 randomly generated 
trees is considered (Fig. 7). All the trees (with their full branch structures) are 
assumed to be of the same species as above, arranged within an approximately 
circular area. The dielectric properties are as previously described. The FDTD 
computational domain has dimensions 29 m x 29 m x 9 m and is discretized into 
1.9 billion cells. Parallelized simulations were performed on a Cray XC30 system 
after partitioning the scene into 256 sub-domains. Each simulation run at 1 set of 


6 











incidence angles required 413 CPU hours. The average response of the forest stand 
at 0 i = 75° is shown in Fig. 8. Imaging of the scene at each elevation incidence 

angle is demonstrated by re-focusing the backscattered fields captured over a 50° 
angular aperture (from (/>. = -25° to 25°, in 2.5° intervals). As seen in 

Fig. 9 (0 i = 75°), the locations of some of the tree trunks can be inferred from the 

co-polarized images. The images, however, are also populated with other clutter 
effects such as multipath interactions (i.e., coupling among trees) and speckle. The 
results at 0 t = 45° and 0 t = 60° are similar to those in Fig. 9. 



Fig. 7 Forest stand with 36 randomly generated trees. Average dimensions: tree height = 
7.3 m; base diameter = 39 cm; canopy diameter = 3.2 m. 
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(a) 



(b) 

Fig. 8 Average response of scene in Fig. 7 
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(a) 



(b) 

Fig. 9 Co-polarized images of scene in Fig. 7: a) vv; and b) hh 
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4. Conclusions 


Full-wave simulations of scattering from trees are considered in this study. 
Geometrically realistic representations of tree structures are first created using an 
open-source random tree generator. Mesh processing and voxelization of the 
structures are then carried out to generate grid models that are amendable to 
numerical analysis. Finally, far-field electromagnetic scattering characterization is 
performed using a parallelized FDTD solver. To demonstrate the practicality of the 
analysis framework, the responses of a quaking aspen tree are investigated as a 
function of frequency, polarization, signal incidence angle, and structural fidelity. 
Subsequently, large-scale simulations are undertaken to examine the responses of 
a forest stand in both the frequency and imaging domains. 
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